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ABSTRACT 

We report the use of the Austraha Telescope Compact Array (ATCA) to conduct 
polarimetric observations of the sky at 5 GHz. The ATCA is normally operated as an 
interferometer array, but these observations were conducted in a split array mode in 
which the antenna elements were used as single-dishes with their beams staggered to 
simultaneously cover a wide area of sky with a resolution of 10/. The linearly polarized 
sky radiation was fully characterized from measurements, made over a range of parallac- 
tic angles, of the cross correlated signals from the orthogonal linear feeds. We describe 
the technique and present a polarimetric image of the Vela supernova remnant made as 
a test of the method. The development of the techniques was motivated by the need 
for wide-field imaging of the foreground contamination of the polarized component of 
the cosmic microwave background signal. 

Subject headings: methods: observational — techniques: polarimetric 



1. Introduction 

Measurements of the polarized component of the cosmic microwave background (CMBP) are 
keenly sought to provide an independent measurement of cosmological parameters and verify the 
precision of certain assumptions that form the basis of cosmological models. Measurements of the 
CMBP can also give greater insights into both the re-ionization and inflation processes than can 
be gained from temperature data alone (e.g. Zaldarriaga, Spergel &: Seljak 1997; Kinney 1999). 
These measurements are best made at frequencies above 30 GHz and are technically difficult. At any 
frequency the measurements are expected to suffer from foreground contaminants such as Galactic 
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synchrotron emission and dust (e.g. Tegmark et al. 2000), with the synchrotron emission expected 
to dominate up to 100 GHz. Observations of Galactic synchrotron emission at lower frequencies 
will provide a measure of its potential to contaminate the CMBP signal. There have been no 
polarization surveys of large areas of sky at frequencies above 2.7 GHz (Duncan et al. 1999) and 
so estimates of contamination at the CMBP observing frequencies currently rely on extrapolations 
from that frequency and from theoretical considerations (Carretti et al. 2002a). Measurements 
of sky polarization at shorter centimetre wavelengths would provide useful constraints on that 
extrapolation. Surveys at higher frequencies have been difficult because of the limited fields of 
view. In general, the need for wide-field imaging has driven the development of interferometers 
with smaller array element sizes and single-dish telescopes with receiver arrays at the focal plane. 
We describe here a novel method of using existing interferometer arrays for wide field imaging. We 
have used the Australia Telescope Compact Array (ATCA), configured as a set of single-dishes, to 
make polarimetric observations at 5 GHz. Single-dish radio-polarimetry has not previously been 
conducted with the ATCA. 

The angular power spectrum of the CMB polarization is predicted to peak around spherical 
harmonic multipoles £ = 1000, corresponding to a ~ 180°/^ ~ 11' angular scale (Zaldarriaga, 
Spcrgcl &: Seljak 1997). At 5 GHz (wavelength of 6 cm), the 22-m ATCA antennas have an 
angular resolution of FWHM ~ 9'. Therefore, the 22-m ATCA antennas, operating in single-dish 
total power mode at 5 GHz, are well suited for making a survey of the sky polarization at the 
multipoles of interest to the planned CMB experiments. 

The sky regions that have the lowest dust emission are the Galactic pole areas. However, 
the minimum Galactic synchrotron emission is observed at lower latitudes (|6| = 40°-50°). In 
particular, a region of interest in the southern sky is at RA = 5'* it l'* Dec = —45° ± 10° (Carretti 
et al. 2002b) which culminates at high elevations at southern mid-latitudes and so is observable 
with the ATCA antennas with low ground spillover. 

We describe the methods used at the ATCA for wide-field radio-polarimetry imaging. The 
Array's standard interferometer mode was used to characterize the instrumental errors and calibrate 
the antenna responses, whereas the wide field imaging was done by using the array in single-dish 
mode with the pointing of the array elements offset to instantaneously cover a wide sky area. The 
method has been demonstrated using the ATCA to image the Vela supernova remnant. 

2. Method 

The ATCA is a six-element aperture synthesis radiotelescope that operates at 20, 13, 6, 3 and 
1.2 cm and 3.5 mm (Prater, Brooks & WhiteoaJi 1992). Each element is a 22-m antenna with 
(shaped) Cassegrain optics. Each has three receiver packages that operate over pairs of wavelength 
bands. Receiver tuning and signal sampling are synchronized across the array by a signal (the local 
oscillator) distributed from the array centre. A range of signal bandwidths can be selected, and 
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each antenna is equipped to provide two observing bands within the tuning range of the receiver 
package in use. In this work we have used the 6/3 cm package tuned to two 128 MHz bands centred 
at 4800 and 4928 MHz. In each band, the signal for each linear polarization is sampled (2-bit) at 
the Nyquist rate, and transmitted on optical fibres to the digital correlator. 

In normal operation the correlator computes the cross-correlation of signal pairs from the 
array, the two polarizations from each of the six elements yielding 60 cross correlations. These 
are calibrated and integrated to form the complex visibilities used in aperture synthesis image 
formation. In addition, autocorrelations and polar cross-correlations are computed for each antenna, 
and in normal operation are used for calibration. We made use of these correlation products for 
the single-dish observations described here. 



2.1. Single-dish polarimetry 

The ATCA antennas are on altitude-azimuth mounts and have two orthogonal linearly polar- 
ized feeds X and Y inclined 45° to the vertical. The radiation field incident on the telescope feed has 
contributions e from the sky, and sources in the immediate environment such as the atmosphere, 
the antenna and the ground. The stray radiation enters the feed both directly and via a number of 
reflections off the telescope optics and support structures. Although initially unpolarized, response 
to this radiation in the X and Y channels may be correlated, resulting in a spurious polarized 
signal. Let s denote the stray radiation and consider it to be partially polarized. The telescope 
optics, even for radiation entering along the intended path, is imperfect so there is some leakage 
of each polarized wave into the other channel. In this section we follow the treatment of radio 
polarimetry by Hamaker, Bregman & Sault (1996), hereafter HBS96. We describe the leakage 
with the terms dx and dy and write a leakage matrix D and input vectors e and s as 



D = 

The signal entering the receiver is then 

D{e + s) = 





dy 1 J \ Gy -\- Sy 



The receiver adds noise {nx, Uy) to each channel and amplifies by gain factors {qx, gy)- Writing 




the input vector v to the correlator is 

V = G[D{e + s)+n] (1) 

^ f 9x \ / ex + Sx + Ux + dx{ey + Sy) \ 
I 9y j \ ey + Sy + Uy - dy{ex + Sx) j 
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The correlator produces the coherency vector which, in the case of single-dish observations, is the 
time-averaged outer product (see HBS96) of the input signal with its complex conjugate 



V = 



/ XX\ 

XY 
YX 



VxV*y 



(2) 



Then (adopting a convention where A^q = apU*) 

XX = g^gliexcl + SxS* + rixnl + dxicyel + Sysl) 

-\- d^ic-xCy + SxSy) 
+ dxdxi^y^y + SySy)] 

~ Gxx\Exx "I" Sxx "I" -^XX "I" dx{Eyx + Syx) 
+ dl{Exy + Sxy) 
+ dxd%{Eyy + Syy)] 

where the <> symbols RIG oniittGci but understood, and terms such as CxsZ, and Sx 

11^, are dropped 

because their factors arc uncorrelated and the time-averaged products vanish. Similarly 

XY = Gxy\Exy ~t~ Sxy ~t~ dxi^Eyy -\- Syy) 
— dy{Exx + Sxx) 
+ dxdy{Eyx + Syx)] 

Note that the receiver noise components Ux and Uy are uncorrelated and so the term < Nxy > 
vanishes and does not appear in the expression for XY. 

For the ATCA antennas, the leakage terms are stable and easily measured with the Array 
configured for interferometry. Moreover, they are small {d< 0.03) and so second order terms {dxd*, 
dxdy, etc.) can safely be neglected for this experiment. The components of V are then 



XX — Gxx[ExX + Sxx + ^XX + dx{Eyx + Syx) + d'^{Exy + Sxy)] 

YY = Gyy[Eyy + Syy + Nyy + dy{Exy + Sxy) + dy{Eyx + Syx)] 

XY = Gxy[Exy -\- Sxy "I" dx{Eyy -\- Syy) dy{Exx "I" "Sjij;)] 

YX = XY* 



(3) 
(4) 
(5) 
(6) 



The quantities of interest arc the Stokes characterisation of the astronomical signal and can 
be related to the terms in {ex,ey) in the expressions above as (see HBS96) 

I + Q' 



Ext. — 
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U' + iV 
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(7) 
(8) 



U' -iV 

Eyx = (9) 

Eyy = (10) 

We seek to measure the linearly polarized sky emission which is characterized by Q and U. In 
Equations 7-10 we use the primed symbols Q' and U' to emphasize that these are the linear Stokes 

parameters in the frame of the antenna feeds, not the required Q and U defined relative to the 
cardinal direction on the sky which must be determined by rotation through the parallactic angle 
tp later in the analysis. Expressed as equivalent temperatures, the typical size of the terms in 
Equations 3-6 for the parameters of our observations are 
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It can be seen that U' can be determined from XY , the component containing terms in Exy^ 
whereas the determination of Q' involves the difference XX — YY. Thus any fluctuations or errors 
in the determination of the receiver gain will, when multiplied by the large terms in Equations 3 and 
4, produce large uncertainties in the value of Q'. For this reason single-dish radio polarimetry is best 
performed with feeds receptive to opposite hands of circular polarization so that all the uncertainties 
arising from differencing the parallel handed components of the coherency vector {RR, LL) flow 
into the circular Stokes V component which is often of less astrophysical interest than the linear 
polarization. Radio interferometry is protected from this because both the stray radiation s and 
the receiver noise n entering the two antennas of each interferometer are uncorrelated. Thus the 
ATCA can successfully use linear feeds for polarimetric interferometry. For this experiment, the 
uncertainties in Q' are intolerable and so we must determine the full linear polarization state from 
measurements of U' made at several parallactic angles. This process is described in section 2.4. 

Further practical difficulties arise from the stray radiation which is partially polarized and is a 
strong function of antenna elevation, and also has some azimuthal dependency. These dependencies 
have proven impossible to model and so all our observations have been conducted as "drift scans" , 
in which the antennas are held at a fixed azimuth and elevation, and the measurements are recorded 
as the sky drifts past at the sidereal rate. 

The analysis above does not include the possibility of further additive components of stray 
signal in the polarized outputs. Examples are the components arising from coupling of receiver 
noise between the two polarizations in the orthomode transducerr, and from common mode noise 
from the use of common Local Oscillators in the frequency conversion stages of the receiver. These 
are expected to be constant, independent of antenna pointing, but are difficult to model. Con- 
sequently absolute measurements of the polarized emission are not attempted. All observations 
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are differential, with images of the polarized emission being presented after the subtraction of a 
baseline from each scan. 



2.2. ATCA amplitude calibration 

At the digitisation stage prior to correlation, the receiver output signals are normalized. At 
the ATCA, the correlator measures correlation coefficients (which take values in the range [-1,1]) 
and to recover the absolute scaling of the coherence measurements a calibration system continously 
measures and records the system temperature Tgyg of each polarization channel against a noise 
signal of known temperature T^ai which is injected at 45° to the X,Y feeds. Typically T^ai — 2K 
referred to the face of the feed horn. The calibration noise signal is switched and a synchronous 
demodulator measures the signal power p and q during the off and on phases respectively. Let a 
be the measurement error on the normalized correlator outputs 



a = 



(11) 



where r and B are the integration time and signal bandwidth, respectively. A factor G can then 
be used to scale the normalised correlator outputs 



GV 



v 



(12) 



q-p 

The measurement error V has contributions from the error AG in the determination of G and from 
the noise a on the normalized correlator outputs 



AG 
~G' 



Aiq-p) 



(13) 



The quantities p and q are independent, being formed by integrating the signal over distinct time 
intervals r/2. Let / = T^al/Tgys- Then q = p{l + f) and 

Ap = V2ap 

Aq = V2aq = V2ap{l + /) 
A{q-p) = ^{Aqf + (Ap)2 



A{q-p) 



2ap ^l + / + p/2 
2a^l + f + P/2 



(14) 



q-p f 

Since at the ATCA / ~ 0.05, the errors in the Vi are dominated by the error in the determination 
of G 



AVj 
Vi 
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/ 



~ 40cr 



when / = 0.05 



(15) 
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Thus, the application of the onhne Tgys cahbration increases the noise in V by a factor of about 
2// and for this experiment severely limits the sensitivity of the measurements. In addition, 
since the calibration signal is injected at 45° to the X^Y feeds, it appears as an additional strong 
polarized component in Equation 5 which, when multiplied by uncertainties and fluctuations in 
Gxy could overwhelm the polarized signal from the sky. Our observations were made with the 
online application of Tcai scaling disabled. The scaling factor G was determined independently 
from separate integrations that bracketted the observations of the sky. 

Note that the difference in the instrumental path lengths for the X and Y polarizations is 
expected to be non-zero, and antenna specific but only slowly changing. Failing to account for it 
rotates the signal between Stokes U and V. We measured this path-length (phase) difference for 
each antenna by measuring the phase of the calibration signal that is injected at 45° to the X,Y 
feeds. The phase correction was then applied to the measured XY values. 



2.3. Observational procedure 

To measure the polarized emission from a Aa x /S.5 rectangular patch of sky centred at (a, 5) 
we performed a series of drift scans of duration At = Aa/ cos 5, with the six antennas pointed at 
the same hour angle and each offset in declination from the next by half the width of the primary 
beam 9. In successive scans the declination was incremented by 30 and the hour angle changed (by 
approximately At) to scan the same range of right ascensions. The entire patch was covered after 

Us = A6/39 

scans, and then repeated for the required total integration time. In this way each point in the 
surveyed area was measured at intervals of approximately 

AaAS 
3d cos 5 



2.4. Determination of Q and U 

As described in section 2.1, the ATCA is able to provide only one correlated output of the linear 
polarization Stokes parameters, U' . Thus Q and U (components defined in sky coordinates) must 
be constructed starting from this single quantity U' by performing at least two scans at parallactic 
angles i/j ideally differing by ±45°. 

The general case consists of N observations of the same pixel performed at several parallactic 
angles. The two Stokes parameters are estimated through a least square approach. Let Ul be the 
z*'* observation of U'. Its expression in terms of the actual Q and U values of the sky in the standard 
reference frame is 

ill = -Q sin(2V'i) + U cos(2 V-i)- (16) 
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where tpi is the parallactic angle during the observation. (In Equation 17, ipi = pi — (f) where Pi is the 
parallactic angle and (j) is the fixed angle between the feeds and the vertical-horizontal directions. 
At the ATCA cp = 45°.) Minimizing the sum of the square differences 

i 

with respect to both Q and U, we obtain as best estimates 



Q = 



E,t// sm{2^|Ji)ZiCosH2^Pi}-J:iU^ cos(2 sm(2 V^,) cos(2 

(Ei sin(2 iji) cos(2 Vi))' - E^ sm2(2 ^Ji) ^ . cos2(2 Vi) 



(18) 



^ ^ EA sm(2^,)E.sin(2^,)cos(2V.,) -E.^/ cos(2 V;,) sin2(2V>, 
(Ei sin(2 Vi) cos(2 V'i))2 - Ei sin2(2 Vi) Ei cos2(2 



(19) 



whose errors, assuming all the samplings have the same sensitivity ai, are 



^2 /„2 



sin2(2 ) [Ei cos2(2 V-i)] ' - E» cos2(2 ^-i) [E^ sin(2 V.) cos(2 



(EiSin(2V',)cos(2V'i))' - E^sin2(2V'^)E^cos2(2V^ 



(20) 



• Ei sin2(2 Vi) Ei sin(2 Vi) cos(2 Vi)]' + Ei cos2(2 V^,) [E, sin2(2 ^,)] ' 



(Ei sin(2 V'i) cos(2 Vi))' - Ei sin2(2 i,,) E, cos2(2 



(21) 



giving the fractional errors with respect to the sensitivity of a single observation. 

These errors depend on the set of parallactic angles at which the measurements occur. The 
efficiency of the reconstruction thus depends on the scanning strategy adopted, which must be 
optimized to minimize the error on both Q and U . Considering that the method provides two 
Stokes parameters from measurements of only one, the ideal sensitivity for N observations is 



,th 



'.V 



(22) 



Its comparison with Equations (20) and (21) allows a quantitative estimate of how well Q and U 
are measured. In addition, to evaluate the overall effect on the {Q, U) pair, we introduce the mean 
error 



2 a? • 



(23) 
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For the simplest case with N=2 observations, the best result is obtained when the difference 
V'l — V'o = i 45°, as expected. 

The case N=3 is less trivial. Keeping in mind the ideal value {alai = 0.82 for N = 3), the case 
is clear looking at Figure 1, where V'o is kept constant, while the other two angles -01 c^nd 02 are 
allowed to vary in the 0°-180° range. The best sensitivity occurs when the three angles are evenly 
separated in a 90° interval, rather than differing by 45°. The case with ±45° separations favours 
one of the two Stokes parameters, resulting in a non-optimal combined sensitivity a. In general, the 
optimum result from N observations is obtained with the 90° range sampled with equally spaced 
parallactic angles. 

3. Test observations 

To provide some confidence in the method, we observed 3.5 x 3.0 degree region centred on the 
Vela supernova remnant for about 10 hr. This region emits strong polarized emission and has been 
well studied at a range of wavelengths and resolutions. Milnc (1980, 1995) reported polarimetric 
observations of the Vela nebula at 2.7GHz (resolution 8.4 arcminutes), 5.0 GHz (4.4 arcminutes) 
and 8.4 GHz (3 arcminutes). Duncan et al. (1997) measured the polarized emission from the 
southern Galactic plane at 2.4 GHz (resolution 10.4 arcminutes), including the Vela region. 

We observed the area with repeated 20-minute drift scans. Each scan was started with all 
antennas set to a = 08'*27"^ and declinations spaced by 5 arcminutes. Thus each scan sampled a 
0.5 degree declination band, and six scans were required to sample the whole area. The six scans 
were repeated for 10 hr, allowing each point in the surveyed area to be measured four or five times. 
The scan data were reduced to form the image shown in Fig. 2. As described in section 2.2, the 
ATCA flux scale is referred to a switched noise signal, the "on-line Tgys measurement". The total 
intensity image was derived from those measurements. The full recovery of Q and U required each 
point to be measured at different parallactic angles (see section 2.4). In this case the south-west 
corner of the surveyed area was observed five times at parallactic angles of —117, —96, —85, —70 
and —48 degrees. 

Figure 2 shows the resulting image of the Vela supernova remnant. The peak total and po- 
larized intensities are 1.3 K and 0.19 K respectively. The image has a resolution of 12 arcmin. 
Our image compares well with published data and verifies our techniques for surveying an area 
in singic-dish mode with series of drift scans and the complete measurement of Stokes Q and U 
through sampling a range of parallactic angles. In particular there is a good match of the polariza- 
tion position angles between our results and those of Milne (1980, 1995). This indicates successful 
removal of the polarization offsets and background polarization which, as reported by Milne (1980) 
is low relative to the polarization of the nebula itself. 

The high brightness of the Vela region makes our test image less useful for assessing the ultimate 
sensitivity of the method and its ability to measure the weak CMBP foregrounds. The apparent 
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noise in the Q and U images of cfq^u — 3 mK is dominated by the variations in polarized emission 
over the field. However, we have made preliminary observations of a region near a = f)^, 5 = —49° at 
high Galactic latitude, which is expected to have low foreground emission. This is the region chosen 
for a number of CMBP measurements (e.g. BaR-SPOrt, Cortiglioni et al. 2003, BOOMERanG- 
B2K, Masi et al. 2005). We made the high-latitude observations using the same method described 
above for the Vela supernova remnant and used a total bandwidth of ~200 MHz. The results of 
the completed 5 GHz ATCA observations of this region will be the subject of a future report, in 
which we will compare them with equivalent measurements at lower frequencies. Reducing the 
preliminary measurements to determine the sensitivity of a single telescope over a 1-s integration 
we find ctq = 5.8 and au = 5.3 mK s^^^, about a factor of 1.6 greater than expected from an ideal 
noise analysis. We can estimate the sensitivity to Q and U for a survey area of f2 with angular 
resolution AO using nt telescopes as 

where T is the total integration time for the survey. For example a 10-hour observation using all 
six ATCA telescopes should yield aq^u — 0-06 mK for a survey area of 1 deg^ and a resolution 
of 12/. This sensitivity would allow a 3cr detection of the expected 5 GHz signal of ~0.2 mK (see 
Carretti et al. (2005) for measurements of the polarized foreground at 2.3 GHz). 



4. Summary 

We have described a novel method, and its demonstration, of using an existing interferom- 
eter array for wide field imaging. The interferometer mode of observing was used to derive the 
instrumental calibration and a single-dish observing mode with pointing offsets between the array 
elements was used for the wide field observing. We used the one correlated output {U') to construct 
both Q and U by observing all points in the imaged area at several parallactic angles. We recognise 
this as a potential technique for all-sky surveying or monitoring with future radio telescope arrays 
such as the Square Kilimetre Array (SKA) if the array elements are chosen to have limited sky 
coverage. This could relax the constraint that the SKA must be built with wide-field elements and 
allow more conventional parabolic dishes to be used as the SKA elements. 
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Fig. 2. — The Vela supernova remnant at 4.9 GHz. The image has a resolution of 12/ (indicated 
at lower left). The 4.9 GHz data were combined with the published 2.4 GHz image (Duncan et 
al., 1997) to determine the rotation measure across the nebula, and thus the intrinsic polarization 
position angles. Total intensity is shown as contours (0.06, 0.10, 0.17, 0.27, 0.44, 0.72, and 1.16 K), 
polarized intensity and magnetic field direction are shown as line segment lengths and angles, 
respectively. The scale of the line segment lengths is indicated by the segment at lower left which 
is equivalent to 0.19 K. 



-13- 



National Facility by CSIRO. 

REFERENCES 

Carretti, E., et al., 2002a in Astrophysical Polarized Backgrounds, AIP Conf. Proc, 609, 109 

Carretti, E., et al., 2002b, in Experimental Cosmology at Millimeter Wavelengths, AIP Conf. Proc, 
616, 140 

Carretti, E. and McConnell, D. and McClure-GrifHths, N. M. and Bernardi, G. and Cortiglioni, S. 
and Poppi, S., 2005, MNRAS, 360, LIO 

Cortiglioni S., et al., 2003, in Warmbein B., ed., 16th ESA Symposium on European Rocket and 
Balloon Programmes and Related Research, ESA Proc. SP-530, p. 271 

Duncan, A.R., Haynes, R.F., Jones, K.L., Stewart, R.T., 1997, MNRAS, 291, 279 

Duncan, A.R., Reich, P., Reich, W., Fiirst, E., 1999, A&A, 350, 447 

Frater, R.H., Brooks, J.W., Whiteoak, J.B., 1992, Journal of Electrical and Electronics Engineering, 
AustraUa - IE Aust. & IREE Aust., 12, 103 

Hamaker, J. P., Bregman, J.D., Sault, R.J., 1996, A&ASS, 117, 137 

Kinney, W.H., 1999, Phys. Rev. D, 58, id. 123506 

Masi S., et al., 2006, A&A, submitted, astro-ph/0507509 

Milne, D.K., 1980, A&A, 81, 293-301 

Milne, D.K., 1995, MNRAS, 277, 1435-1442 

Tegmark, M., Eisenstein D.J., Hu W., de Oliveira-Costa, A., 2000, ApJ, 530, 133-165 
Zaldarriaga, M., Spergel, D.N., Seljak, U., 1997, ApJ, 488, 1 



This preprint was prepared with the AAS IM^^X macros v5.2. 



